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Abstract 

We study the 50(3) lattice gauge theory in 3+1 dimensions with the adjoint Wil- 
son action modified by a Z2 monopole suppression term and by means of the Pisa 
disorder operator. We find evidence for a finite temperature deconfinement transi- 
tion driven by the condensation of U(l) magnetic charges. A finite-size scaling test 
shows consistency with the critical exponents of the 3D Ising model. 



Lattice SU (N) pure gauge theories in the fundamental representation undergo 
a finite temperature deconfinement phase transition [1,2] signalled by the spon- 
taneous breaking of the global center symmetry Zjy [3,4]. It is an interesting 
question whether this happens also for the theory in the adjoint representa- 
tion SU(N)/7*n [5]. In the latter case the center is trivial and naively there is 
no global symmetry to be broken. Moreover, although both the fundamental 
and the adjoint SU(N) lattice theories have the same naive continuum limit, 
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non-perturbative investigations carried out for Wilson [6,7] as well as for Vil- 
lain discretizations [8] showed a different behavior. For the adjoint cases no 
evidence for a finite temperature phase transition was found, whereas a bulk 
transition separating the strong from the weak coupling regions appeared for 
both kinds of discretizations. 

The bulk transition was explained in terms of a condensation of lattice arti- 
facts - Z 2 monopoles [8] . It was argued that the finite temperature transition 
could be overshadowed by the bulk one and Z 2 monopole suppression terms 
were thus proposed [9]. More recently the Villain mixed fundamental-adjoint 
SU(2) model with a monopole (and vortex) suppression has been reinvesti- 
gated [10,11,12] and first evidence for a deconfinement transition in the Ising 
3D universality class, at least for strong coupling (JV r = 2), was given. More- 
over, in the weak coupling region negative states of the Polyakov loop in the 
adjoint representation were found [10,13] and later linked to the non-trivial 
twist sectors of the theory [14], given that on the hypertorus T 4 the parti- 
tion function of the 5*0(3) theory in the Villain formulation with complete 
Z 2 monopole suppression was shown to be equivalent to that of the SU{2) 
theory in the fundamental representation when summed over all twist sectors 
[15,16,17,18]. 

In a recent paper we have reinvestigated the 5*0(3) theory using the Wilson 
action and a Z 2 monopole suppression term with a "chemical potential" A 
[19,20,21,22]. The phase diagram in the (3a — A plane was studied at zero 
and finite temperature monitoring the tunneling between twist sectors and 
its disappearance for strong enough Z 2 monopole suppression. For A > .85 
we have found a strong indication for the existence of a finite temperature 
deconfinement transition although having restricted the simulation to a fixed 
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(e.g. trivial) twist sector. The proposed phase diagram is redrawn in Fig. 
1. Since no proper order parameter was available a determination of critical 
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Fig. 1. Phase diagram of the modified 50(3) Wilson theory with %2 monopolc 
suppression as seen in the trivial twist sector. The phases are characterized by 
(non)vanishing values of the monopole density M, the average adjoint Polyakov 
loop (La) and the averaged electric twist variable (z) (see [21]). 



exponents was not intended, and there was no answer given to the question 
about the underlying confinement mechanism. 

A disorder parameter related to Abelian monopole condensation in the dual 
superconductivity picture of confinement [23,24,25] has been devised by the 
Pisa group some time ago [26,27,28,29,30]. It is the vacuum expectation value 
of a magnetically charged operator (//) shown to be different from zero in 
the confined phase, thus signaling dual superconductivity, and going to zero 
at the deconfining phase transition. Similar parameters have been constructed 
more recently by Frohlich and Marchetti [31] as well as in the framework of 
the lattice Schrodinger functional [32,33] leading to analogous results. The 
main advantage of these parameters is that they can be applied also to full 
QCD, where the center symmetry is explicitly broken by the fermionic degrees 
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of freedom and - as we shall show - to the adjoint pure gauge theory, where 
center symmetry becomes trivial. In this letter we will use the Pisa disorder 
operator in order to answer the questions raised above for the SO (3) lattice 
gauge theory. In particular we will check the critical exponents and whether 
the dual superconductor scenario applies also to this case. 

We will study the SU{2) adjoint representation Wilson action modified by a 
Z 2 monopole suppression term 

^^(l-^j+^-c), (1) 

where the product a c = n Pedc sign(Tr pUp), taken around elementary 3-cubes 
c, defines the Z 2 magnetic charges. Their density can be introduced as M = 
1 ~~ (l^~I]c "c) normalized such that it tends to one in the strong coupling 
region and to zero in the weak coupling limit, N c denoting the total number 
of elementary 3-cubes. Although a c is constructed in terms of fundamental 
representation plaquettes, it is a natural SO (3) quantity ensuring that the 
action (1) is center-blind in the entire Pa — A plane. The link variables can 
be represented both by 5*0(3) or SU(2) matrices, exploiting the property 
Tr^ = Tr^ — 1 for the Wilson term or picking a random SU (2) representative 
of the 50(3) link to construct the Z 2 monopole contribution. A standard 
Metropolis algorithm has been used to update the links. In [21] we have found 
a strong indication in favor of finite temperature phase transitions with lines 
moving up with the time-like lattice extent N T and running away from the 
bulk transition approximately parallel to the A axis (see Fig. 1). 

The Pisa disorder operator [26,27,28,29,30] was shown to be a reliable order 
parameter for SU{2) and SU{3) gauge theories in the fundamental repre- 
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sentation, with and without dynamical quarks, giving critical exponents in 
agreement with other order parameters. Moreover, it can give important in- 
formations about the mechanism which confines quarks into hadrons. The 
Pisa disorder operator is motivated by the dual superconductor scenario for 
the QCD vacuum [23,24,25], driven by the condensation of U(l) magnetic 
charges. The construction of the operator in the case of the modified 5*0(3) 
theory follows the same line as in the fundamental case, so we will avoid go- 
ing into the details and we will refer to the original papers for further details 
[26,27,28,29,30]. 

The idea is to construct a magnetically charged operator /i which shifts the 
quantum field at a given time slice by a classical external field corresponding 
to a magnetic monopole. The U(l) subgroup of the gauge group which de- 
fines the magnetic charge is selected by an Abelian projection, usually fixed 
by diagonalizing an operator X in the adjoint representation. The disorder 
parameter is defined as 

<^> = «^ 

where Sm^) denotes the Wilson action with the space-time plaquettes U^x, t) 
at a fixed time-slice t modified by an insertion of an external monopole field 

U i4 (x,t) = U i (x,t)$ i (x + ly)U 4 (x + lt)U}(x,t + l)Ul(x,t)) : (3) 

where <&i(x, y) = ile iTab i( x ~ l '^& , with f2 the gauge transformation which diag- 
onalizes an operator X in the adjoint representation. T a denote the generators 
of the Cartan subalgebra and b the discretized transverse field generated at 
the lattice spatial point x by a magnetic monopole sitting at y. We decided to 
work with the completely random Abelian projection (RAP) [29,30] in which 
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we do not diagonalize any operator X: it can be thought as a kind of averaging 
over a continuous infinity of Abelian projections. It should be stressed that 
only the plaquette contribution to the action (1) is modified by the insertion 
of the monopole field and not the chemical potential term. From the definition 
of /i, making use of an iterated change of variables, it can be shown that the 
correlation function D(At) = (fi(y,t + At)/i(y,t)) describes the creation of a 
monopole at (y,t) and its propagation from t to t + At [26,27,28,29,30]. At 
large At, by cluster property, D(At) ~ Aexp(-MAt) + (/i) 2 . A non-vanishing 
(/i) indicates spontaneous breaking of the U{1) magnetic symmetry and hence 
dual superconductivity. In the thermodynamical limit one expects (/i) ^ 
for T < T c , while (/i) = for T > T c if the deconfining phase transition is 
associated with a transition from a dual superconductor to a trivial vacuum. 

At finite temperature there is no way to put a monopole and an antimonopole 
at large distance along the t-axis as it is done at T = 0, since at T ~ T c 
the temporal extent N T a is comparable to the correlation length. Therefore, 
one computes (/i) but with C*-periodic boundary conditions in time direc- 
tion imposed to the numerator in Eq. (2) in order to ensure magnetic charge 
conservation: Ui(x,N T ) = U*(x,0), where U* is the complex conjugate of Ui. 
These boundary conditions have been indicated in Eq. (2) by the index M at 
the integration measure. They change the sign of the term proportional to cr 3 
in the links, creating a dislocation with magnetic charge -1 at the boundary 
which annihilates the positive magnetic charge created by the operator /i . 
As a consequence the magnetic charge is conserved and everything is con- 
sistent. Of course, the denominator in Eq. (2) is computed using standard 
periodic boundary conditions. We used links in the fundamental representa- 
tion, so we implemented exactly the above condition, but an analogous con- 
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dition holds also for link variables defined in the adjoint representation, i.e. 
Ui(x, N T ) = (I 3 + 2T 2 2 )[/j(a;, 0)(I 3 + 2T 2 2 ); charge conjugation is realized in both 
representations through rotations by an angle 7r around the color 2-axis. 



Since (//) is the average of the exponential of a sum over the physical volume, 
it is affected by huge fluctuations which make it difficult to be measured in 
Monte Carlo simulations. A way out is to compute the derivative with respect 
to the coupling parameter (5 = (3a, which contains all the relevant information 



It is given by the difference between the Wilson plaquette action term (IT) 
averaged with the usual measure and the modified plaquette action term (Um) 
averaged with the modified measure {DU) M e~ SM I J(DU) M e- SM . The order 
parameter can be reconstructed from 



Eq. (5) tells us that in order to have (p) ^ in a confined phase, where 
the dual magnetic symmetry becomes broken, p should stay finite for j3 < 
f3 c in the thermodynamical limit (i.e. in the limit of spatial lattice extent 
N s — > oo), while a sharp negative peak for p occuring at f3 c and diverging for 
N s — > oo should signal the phase transition associated with the restoration of 
the dual magnetic symmetry. Above the transition a sufficient condition for 
(p) to vanish would be to have p — > — oo in the thermodynamical limit. In 
[26,27,28,29,30] for the SU(2) case it was argued using perturbation theory 
that p for (3 — > oo reaches negative plateau values linearly scaling with N s . In 
[34] a more detailed numerical analysis has been performed for both SU(2) and 
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\og(p) . 



(4) 
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SU(3) pure gauge theories, showing that, in the case where /i is magnetically 
charged, p diverges (to negative values) linearly in the weak coupling limit, 
and more and more rapidly as T — > T c from above, where it diverges as N^ u , 
thus proving that (p) is strictly zero for every temperature T > T c , as follows 
from Eq. (5). Therefore, when simulating the theory at accessibly large j3 
values and lattice sizes, a good criterion for (p) being exactly zero above the 
transition temperature, is that p keeps diverging at least linearly with N s in 
a wide range of (3 values above the transition. 





(a) A = 0.0 



(b) A = 0.4 





(c) A = 0.7 



(d) A = 0.8 



Fig. 2. Disorder operator p computed at different values of the chemical potential 
A at finite temperature in the asymmetric volume V = 4 x 12 3 . 
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We will now present our numerical results. First of all we studied the bulk 
transition (compare with Fig. 1) for some A-values between (no monopole 
suppression) and 0.8 (partial monopole suppression) by varying (3a- We do 
this for finite temperature (N T = 4) to investigate the interplay between the 
bulk transition and the finite temperature one. Fig. 2 shows a clear dip at the 
location of the bulk transition, in agreement with what was found in previous 
works with other observables. Although the dip in p decreases in magnitude 
with increasing A, i.e. as we are suppressing more and more the Z 2 monopoles, 
the lattice artifacts are still present and overshadow the finite temperature 
transition, assuming the latter exists. Therefore, one cannot yet see neither 
the scaling of the physical transition with the temperature nor the finite-size 
scaling of the operator with the spatial volume, although one expects that on 
large enough lattices the two transitions should decouple also for these values 
of the couplings [21]. Anyway we can conclude that a condensation of U(l) 
magnetic charges takes place below the bulk transition line. 

Let us now turn to the more interesting case of stronger monopole suppression 
(A = 1.0), where the decoupling of the bulk from the finite temperature tran- 
sition occurs already at reasonable volumes. Again we vary the temperature 
through (3 a- We have kept the system fixed in the trivial twist sector so to 
compare with the results for the adjoint Polyakov loop and the distribution 
of the fundamental Polyakov loop variable in [21]. As one can see from Fig. 
3 (a), at fixed N T = 4 the parameter p shows a dip around (3a — 1.0 which 
becomes deeper and deeper with increasing spatial volume, p stays finite in 
the low Pa region as the spatial volume is increased, as can be inferred from 
Fig. 3 (c), while the data for high (3a are consistent with its saturation at neg- 
ative plateau values diverging more than linearly with N s , as Fig. 3 (d) shows. 
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Fig. 3. p computed in the trivial twist sector at finite temperature (N T = 4) for dif- 
ferent values of the spatial volume and for chemical potential A = 1.0 (a). Finite-size 
scaling analysis for p (b). The strong and weak coupling regions are highlighted in 
(c) and (d) respectively. 



This gives us clear evidence that also in this case, like for the lattice theory in 
the fundamental representation, condensation of U(l) magnetic charge takes 
place below the deconfmement phase transition and disappears above T c . The 
vacuum is a dual superconductor in the confined phase and becomes trivial in 
the deconfined one. It shows that the Pisa disorder operator is a meaningful 
order parameter also in a center-blind theory and for fixed twist. Moreover, 
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the position of the dip is well-consistent with the results found in our previ- 
ous work [21]. The statistics is not sufficient in order to determine the critical 
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Fig. 4. p computed in the trivial twist sector at finite temperature (iV r = 6) for 
different values of the spatial volume and chemical potential A = 1.0 (a). Finite-size 
scaling analysis for p (b). 



exponents independently, but as one can see from Fig. 3 (b), using the known 
critical exponent for the 3D Ising model, v = 0.63 and estimating the critical 
adjoint coupling (3 C A = 0.98 the data show a reasonable finite-size scaling be- 
havior with the spatial volume. We have also checked the scaling with respect 
to the continuum limit for varying N T , making some simulations by fixing 
N T = 6 and N s = 16, 20, 24. As one can see from Fig. 4 (a) the dip occurs 
at a larger (3a value than for N T = 4 and becomes again deeper by increasing 
N s , in agreement with a finite temperature phase transition. In this case we 
estimate (3 C A = 1.19. Fig. 4 (b) shows the quality of finite-size scaling assuming 
the Ising model value for the critical index. Our results are in agreement with 
the finite-size scaling observed for the specific heat in the Villain action case 
for N T = 2 [10,11,12]. 

We can conclude that our investigation of the modified lattice 5*0(3) gauge 
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theory with Wilson action and Z 2 monopole suppression using the Pisa dis- 
order operator has confirmed our previous results [21]. Although having re- 
stricted to the fixed trivial twist case we find a clear indication for the ex- 
istence of a finite temperature transition decoupled from the notorious bulk 
transition. The critical behavior at such a transition reasonably agrees with 
the critical exponents of the 3D Ising model. The nature of the Pisa disorder 
parameter we used tells us that the transition is related to a condensation 
of U(l) magnetic charges for T < T c which disappears above the transition, 
proving that the dual superconductor scenario is a good model of confinement 
also for the adjoint theory. Of course, a final answer can be given only after 
all twist sectors will be taken into account simultaneously. This study is cur- 
rently under way. Moreover, it would certainly be interesting to compute also 
the free energy of an (extended) center vortex in order to check how the vortex 
condensation mechanism works in this case. 

We thank A. Di Giacomo, L. Del Debbio, and P. de Forcrand for valuable 
comments and discussions. 
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